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Abstract 

In the experiments on stress-induced phase transitions in SMA strips, several 
interesting instability phenomena have been observed, including a necking-type in- 
stability (associated with the stress drop), a shear- type instability (associated with 
the inclination of the transformation front) and an orientation instability (associ- 
ated with the switch of the inclination angle). In order to shed more light on these 
phenomena, in this paper we conduct an analytical study. We consider the problem 
in a three-dimensional setting, which implies that one needs to study the difficult 
problem of solution bifurcations of high-dimensional nonlinear partial differential 
equations. By using the smallness of the maximum strain, the thickness and width 
of the strip, we use a methodology, which combines series expansions and asymp- 
totic expansions, to derive the asymptotic normal form equations, which can yield 
the leading-order behavior of the original three-dimensional field equations. An im- 
portant feature of the second normal form equation is that it contains a turning 
point for the localization (necking) solution of the first equation. It is the presence 
of such a turning point which causes the inclination of the phase front. The WKB 
method is used to construct the asymptotic solutions, which can capture the shear 
instability and the orientation instability successfully. Our analytical results reveal 
that the inclination of the phase front is a phenomenon of localization-induced buck- 
ling (or phase-transition-induced buckling as the localization is caused by the phase 
transition). Due to the similarities between the development of the Luders band in a 
mild steel and the stress-induced transformations in a SMA, the present results give 
a strong analytical evidence that the former is also caused by macroscopic effects 
instead of microscopic effects. Our analytical results also reveal more explicitly the 
important roles played by the geometrical parameters. 
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1 Introduction 



Shape memory alloys (SMAs), due to their two important characteris- 
tics, shape memory effect and pseudoelasticity, have broad applications (see 
Duerig et al 1990, Otsuka and Wayman 1998). To understand the behaviour 
of this type of materials, systematic experiments have been carried out on 
uniaxial tension of superelastic NiTi (a kind of SMAs) wires, strips and tubes 
(see Shaw & Kyriakides 1995, 1997, 1998, Sun et al 2000, Tse & sun 2000, 
Feng & Sun 2006, Chang et al 2006). An experimental vedio on the ten- 
sion of a NiTi strip by Q. P. Sun's group can be found in the website: 
www.me.ust.hk/meqpsun/vedio/tension-strip.htm. Among many important ob- 
servations in these experiments, some key features are the various instability 
phenomena associated with stress-induced phase transitions. For example, for 
the stress-induced phase transitions in a strip during a loading process, at least 
three instability phenomena were observed (see Shaw & Kyriakides 1998, Sun 
et al 2000, Tse & Sun 2000): (i) a stress drop after the nucleation of the marten- 
site phase and associated with it there is a formation of two phase fronts, 
which manifest like a neck (a necking- type instability); (ii) the phase front 
inclines an angle with the strip axis (a shear- type instability); (hi) the front 
can switch the inclination to an opposite angle. Finite element simulations 
have been carried out (see Shaw & Kyriakides 1998, Shaw 2000) to capture 
the main features observed in experiments. The numerical results in these two 
papers revealed some important information about the stress-induced phase 
transitions in strips. For example, it was found that the evolution of phase 
transition events is strongly influenced by overall geometric (structural) ef- 
fects. The results of these two papers strongly suggest that continuum level 
events remain dominant players in the SMAs considered by them. 

Motivated by the experimental and numerical results by others mentioned 
above, in this paper we shall study instability phenomena during the phase 
transitions in a strip analytically. We model this problem in a continuum three- 
dimensional setting with a no-convex strain energy function, in view of the 
results of Shaw & Kyriakides (1998) (cf. the last sentence of the paragraph 
above). The difference between the stress-strain relation used in this paper and 
Shaw & Kyriakides (1998) is that the former is a cubic nonlinear curve while 
the latter is a trilinear one. The nonlinearity could play certain role (see, Fig. 8 
of Shaw 2000). Since in the experiments, the maximum strain is less than 8%, 
keeping the nonlinearity up to the third order (see (2.3)) is accurate enough, 
at least not worse than a trilinear approximation. We also point out that here 
the intention is to study macroscopic instability phenomena in the loading 
process only and no attempt is made to consider the microscopic effects. 
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Analytical results, if achievable, have a number of advantages. One is that 
there is no need to introduce an artificial imperfection to capture the post- 
bifurcation mode. Secondly, mathematically, an instability is caused by the 
fact that there are multiple solutions, and analytical results can shed light 
on how this situation arises and help to understand the mechanism. Thirdly, 
from the analytical results one can see more clearly the roles played by various 
parameters (for the present problem, in particular, the geometric parameters, 
e.g. the thickness and the width). Indeed, the analytical results obtained in 
this paper reveal more explicitly the important role of the thickness of the 
strip and show that the width influences the instability phenomena through 
the thickness-width ratio rather its magnitude. 

As pointed out in Shaw & Kyriakides (1998), in the macroscopic scale there 
are many similarities between stress-induced transformations in a SMA and 
the development of Luders bands in a mild steel. In the literature, there are 
different views whether the Luders band is caused by microscopic effects or 
macroscopic effects (see Estrin & Kubin 1995). Here, we have shown that the 
inclination of the transformation front is a phenomenon of localization-induced 
buckling. This offers a strong analytical evidence that this phenomenon in a 
SMA, and plausibly the phenomenon of Luders bands (due to the similarities), 
is due to macroscopic effects. 

Since we formulate this problem in a three-dimensional setting with a nonlin- 
ear constitutive relation, the governing field equations are three coupled non- 
linear partial differential equations (PDEs). It is extremely difficult to deduce 
the post-bifurcation solutions of nonlinear PDEs analytically. Fortunately, for 
the present problem, several quantities are small, e.g., the thickness, the width, 
and the maximum strain, which then permit us to use a methodology of cou- 
pled series and asymptotic expansions to deduce asymptotic solutions. This 
methodology was first introduced to study nonlinear waves in solids (see Dai 
& Huo 2002, Dai & Fan 2004). Recently, it has been successfully used to study 
various instability phenomena in solids (see Dai & Cai 2006, Cai & Dai 2006, 
Dai et al 2008, Dai & Wang 2008). However, all those problems studied before 
are essentially two-dimensional. Here, for the first time this methodology is 
used to study a three-dimensional problem. 

The remaining of this paper is arranged as follows. In section 2, we give the 
general three-dimensional field equations for a plate. Then, in section 3, we 
non-dimensionalize the three-dimensional governing equations to identify the 
key small variables and small parameters for a thin plate. And then, by using 
the smallness of these variables and parameters and a methodology of coupled 
series and asymptotic expansions, in section 4 we derive the asymptotic two- 
dimensional equations. By considering the smallness of the width further and 
using a similar methodology in section 4, we obtain the quasi one-dimensional 
asymptotic normal form equations for a thin strip in section 5. These normal 
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form equations are then solved analytically in section 6 for an infinitely long 
strip and the solutions obtained seem to be able to describe many features 
observed in experiments. Finally, some conclusions are drawn. 



2 Three-Dimensional Field Equations 



We consider the deformation of a three-dimensional plate composed of 
a hyperelastic material. Its thickness is 2a and its width is 2b. We use the 
Cartesian coordinates (x, y, z) (equivalently Xi) and (X, Y, Z) (equivalently 
Xi) to represent a material point in the current and reference configurations, 
respectively. The geometry of the object of study is shown in Figure 1. 




Fig. 1. The geometry of the object of study 

To ensure phase transitions can take place, we suppose that the strain energy 
function $, which is a function of the invariants of the left Cauchy-Green 
strain tensor for a homogeneous isotropic hyperelastic material, is non-convex 
such that there is a local maximum and a local minimum in the uniaxial 
stress-strain curve under a homogenous constant strain state. The first Piola- 
Kirchhoff stress tensor S is given by 



S = ^-, (2.1) 
where F is the deformation gradient and the components of (F — I) are 



Ou ■ 

UtJ= dX~> 6.7 = 1,2,3), (2.2) 

where Ui are the components of the displacement vector and Xj are coordinates 
in the undeformed configuration. If the strains are small, it is possible to 
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expand the first Piola-Kirchhoff stress components in term of the strains up 
to any order. Due to the complexity of the problem, we shall consider the 
material nonlinearity up to the third order. The formula containing terms up 
to the third-order material nonlinearity has been provided in Fu and Ogden 
(1999) as 



^ji — a jilk U k,l + 2 a jHknm U k,l u m,n + Q a %iknmqp U k,l u m,nUp jq , (2.3) 

where ah ilk , a^n knm an d a %iknmq P are incremental elastic moduli, which can be 
calculated once a specific form of the strain energy function is given. Their 
expressions can be found in Appendix A, where it is also shown that a^ ilk 
has 4 non-zero members and only two are independent, Oj ilknm has 9 non- 
zero members and only three are independent and Oj Uknmgp has 22 non-zero 
members and only four are independent. For the convenience of the sequel 
analysis, we write out the index 3 explicitly. As a result, from Eqs. (12.31) we 
have 



^ji - a )i33 M 3,3 + al jia3 U 3,a + «jj 3 a M a,3 + a )if3a U a,l3 

/ 2 2 / 2 2 2 \ 

+ 2 ^j^S^S + 2 ( a jia333 M 3,a + «-jj 3a3 3 M a,3 + a ji/3a33 U a ^) U 3i3 

2 2 2 

+ a jia3f33 U 3,aU 3 ,f3 + 0, ji3a3j3 U^ 3 Up^ 3 + a ji(3aay U a,l3Uy,a 

+ 2 ( a2 ji3/3a3 U 3,aUp,3 + 0>%/3ar t 3 U a,pU3,'y + ^^U^pU^^)) 

+ g {3{ a ji@a3333 u a!,/3 + a ficc33333' u 3,a + a %3a3333 u ot,3) u 3,3 
+ 3 ( a jia3/3333 U 3,a«3,/3 + a ^3a3/333 U «,3«/3,3 + a ji/3^733 M a,/3 M 7^) M 3,3 
+ §(. a )i3i3a333 U 3,a.Up > 3 + a ji/3a3 7 33 U a,/3 W 7,3 + a jj/3a<7333 U a,/3M3, 7 )«3,3 
+ a jia3/33j3 U 3,aU3,pU 3a + (a^^^Ua^Up^ + 2,a? jia3 p 33l U 3 ^ a U 3 ^ 

+ 3a fia33/33 7 ' U 3,a«A3)'U7,3 + {^]iaSP3^ U 3,aU3,P + ^%Za3^ U a,3Up,3 
+ ^ >a ]iZaP3a 1 u a,3U3,p + 3a 3 ji3rf3a(T1 U Tj3 U at(3 + 3a^ iT30aa7 U 3)T U a> p)u 7>a 
+ a ji333333 M 3,3 + tfirXPav^x^Ua^U^) , (a, (3, 7, (7, T, A = 1, 2). 

(2.4) 

The equations of equilibrium are 



E iw = 0. (2.5) 
Substituting f!2.4p into (12.51) . we obtain 

^333^3,33 + (bi3 a 3 + bi a33 )u 3)3a + &i/3 a3 M3,a/3 
+ fri33a^a,33 + (^i/33a + bi3/3a)u a ,3/3 + h^PaUa,^ = 0, (2.6) 
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where 



7 1_L 2 _l_ 2 _l_ 2 _l_ 2 

Oijkl — ttjilk ~r a jilk33 U 3,3 i a jilkf3a U a,f3 ~+~ a jilk3a U a,3 > a jilka3 U 3,a 

1 g 2 / 3 3 3 \ 

+ T^j^SSS^S.S + ( a jilkf3a33 U a,P + ^ii/caSSa^a + a ji/fc3a33' U a,3 J w 3,3 

+ lj a %ikf3ao~i U ot,l3 u ~j,cT + ( a jiifc«3<T7 M 3,o + a %ik3aa-y U a,3)u liCT 

+ a3 jilk/333a U a,3U3,l3 + ^{^jilkpZa^^a^P + 3 jilk3/33a U a,3Uf3,3) ■ (2-7) 

It should be noted that, for an initially isotropic material these terms in the 
right hand side of (12.71) will vanish when the number of index 3 in their all 
subscripts is odd. 

Eqs. (12. 6p are the governing equations for the three unknowns Ui(i = 1, 2, 3). 
In order to investigate the instability phenomena, one needs to study the solu- 
tion bifurcations of the three-dimensional nonlinear partial differential equa- 
tions (PDE's) (12. 6p with the traction free conditions on the top/bottom sur- 
faces and two side surfaces and under proper end conditions. Mathematically, 
this is a very challenging problem, since there is no available a general method 
for studying the bifurcations of three-dimensional nonlinear PDE's. Here, we 
shall use a novel approach involving coupled series and asymptotic expansions 
to derive the asymptotic normal form equations in order to carry out the anal- 
ysis. For that purpose, we first non-dimensionalize the governing equations to 
identify the key small variables and small parameters. 



3 Non-dimensional Equations 

Suppose that the loads acting on the boundaries of the plate are symmetrical 
about the mid plane and therefore the deformation is also symmetrical about 
this plane. Then, we have 

u a {X p , -X 3 ) = u a (Xp,X 3 ), u 3 (X p , -X 3 ) = -u 3 (X p ,X 3 ). (3.1) 
Based on Eq. (13.11) . we introduce a transformation 

u 3 = X 3 w, s = Xj. (3.2) 
The dimensionless quantities are defined through the following scalings: 

s = l 2 s, X a = lx a , u a = hu a , w = —w, (3.3) 
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where / is the length of the plate and h is a characteristic displacement in the 
mid-plane. From (13. 2p and (13.31) . we obtain 



„dw du r 
u 3 ,3 = e{w + 2s— ), u a> p = e- 



OS OX/3 

U3, a = eV§—^, u a>3 = 2e\/S—^, (3.4) 
ox a OS 

where e = h/l is a small parameter (i.e., we are considering a weak nonlinear- 
ity). The second-order derivatives can be treated similarly. 

Substituting (JS3D and ([33]) into Q22D and (T2SD, we obtain 



a PT-ya U a,l3'y + 2a 3 T 3 a U a,s + ( a 3ra3 + a ar 33 ) W ,< 



+ s(4a3r 3a M a,^ + 2 ( a 3ra3 + a l 



ar 33 ) W ,as 



+ e { a2 f3T- 1 ar l S, u a,l3^ u ^,ri + 2ci3 r 3 m? gWa,sW§,» ) + Op^^Ua^W + 2al Ta333 U ajS W 

HalraSvH + a ar33^K,^,a + (<4ra333 + a lr 3333 )^W , a 
+ S (^( a 3raf3rj3 + a 'pT 3 ar] 3 ) U a,l3sUr 1 ,s + ^a^g^liQ^glt,,^ 

+ 4a 3™333( M «,^ w + V,,) + 2a 2 aT03rj3 u v , s w ta p + 2a 2 p Tia33 u a:J3j w tS 
-\-12a 3T33a3 u a)S w jS + 2(a 3TQ , 3 33 + 0^3333 + 3a 3r333a )w ta w iS 

+ 2 ( a 3r3 7 /3a + a ?r33/3a K,/? W ,7S + 2 Kr3333 + a l T1 333) WW ,is 
+ 2 (a 3T /3a3r, + ^rSaSr,) U a ,p a W lT , + a 2 QT , a pW :V ) 

+ As 2 {2a\ Ta333 U a ^ s W tS + (a„ T 3333 + a 3 Va333) U 'a,s U ',.s 
+ 2a 3r3 3 Q , 3 Ua,sti',ss + °3T333a W ,a' u; ,ss)) 

+ e 2 (tf 1 ) = (r = l,2), (3.5) 



2 ( a 33/3a + a J 3 3a>a,/? S + °a3/33™ a/3 + 603333^ + 4s<4 333 W sfl 

^2 /i2 ^2 2 \ 

+ e { 2a i33'yari3 U a,f3'y u ri,s + ^a^s^Ma^M^s + 2(a 33/ g Qr) £ + a a3 ps rl (:)'U>a,f3s u £,r) 

+2(a 33( g a 33 + aQ,3/33 3 3)w ai/ g s W + 2(0^3333^ + 0/9333^3)^,5^^ 

+^3/33^,*?™ <*/3 + ^3/3333^,^/3 + 6a 3333/3a U a,/3™ s + 6^33333 

Jra 'p3~/a3r l U a,f3'y w ,n + 2a^, 3 333 T? 'U aiS W i?? + 20^,3333^10^^^ 

/ 2 / 2 2 \ 2 

+ 2s(4a Q , 3 33 /33 'U Q , iSS -U ( g jS + 2(a 33/3a 33 + a Q ,3 ( g333)Ma,/3sW jS + 00.3/3333 W a /3«; )S 

+ a 333333( 2 ^,ss + QW 2 S ) + A(X 2 a33W3 U^ s W, as + 2a^33 3/ 3 Q M ai /3W >ss 
+2 a a3333/3 M a,ssW ; ,/3 + ( a a 3 333/3 + a a3333/3) U ',/3 u ',as) 

+ 8s 2 a 333333 w tS w tSS ) 

+ e 2 (H 2 ) = 0, (3.6) 
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where and thereafter the tilde over non-dimensional variables has been dropped 
for convenience. The lengthy expressions for Hi(i = 1,2, • • •) are omitted al- 
though they are needed for the calculations (interested readers can contact 
the corresponding author for their expressions). We consider the case that the 
top and bottom surfaces of the plate are traction- free. By using (|2.4p . we have 



3a W ,a 



^ a 3ra3 U a,s + a 3 T ; 

+ e03ra/3 7 3 M /3,^7,s + 4 °3™333 M «,^ + 2a 3™/33 7 M /3,« W i 7 + 2a 3r3a33 WW ,/3 

+ 2sw !S (2al Ta333 u a)S + a\ T3a33 w <a )) 

+ e 2 (H 3 )\ s=Ul =0 (r = l,2), (3.7) 

«33a/3«/3,a + 4j33( W + 2sW ,s) 

C j 2 o 2 2 2 

+ 2^ a 33a/3^/6 u l3,aU8,-y + 2 a 33al333 U i3,a w + ^SSSS" 10 

+ s(4al 3a . m u a , s Uf3, s + 4a 2 33al333 u^ a w, s + 4a 2 333333 ww, s 

+ 4a 33a33/3 W a, S ™ /3 + a 2 333a3f3 W yQ W )t3 ) + A S 2 ^333333^) 

+ e 2 (# 4 )U„ 1 =0, (3.8) 

where v\ = a 2 /I 2 . Please note that due to symmetry, the boundary conditions 
at the bottom surface are automatically satisfied. 

Eqs. (13.51) and (13.61) provide the governing equations for three unknowns u a 
and w and the boundary conditions are (13.71) and (13. 8ft . However, they still 
comprise a formidable system of nonlinear PDE's to be analyzed directly. To 
go further, we assume that the plate is thin. Then v\ is a small parameter and 
< s < Vi is a small variable. It is clear that the unknowns are functions of 
the spatial variables X\ and x^, the small variable s and two small parameters 
Vi and e. Next, we shall use the smallness of the variable s and two parameters 
Vi and e to derive the asymptotic two-dimensional equations. 

Remark: Since the current methodology to deduce the one-dimensional asymp- 
totic normal equations from the three-dimensional nonlinear field equations 
has not been done before, in the next two sections we shall provide some 
detailed derivations. 



4 Two-dimensional Asymptotic Equations 

As discussed in the previous section, we can write 

u a = u a (x /3 ,s;e,v 1 ), w = w(x p , s; e, v l ). (4.1) 
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As the variable s is small, as long as we assume that the unknowns are suffi- 
ciently smooth in s, we can take the series expansions in s for the unknowns, 
i.e., 

u a = Uoaixp; e, v x ) + sU la {xp] e, Vx) + s 2 U 2a (x p ; e,v x )-\ , 

w = W (xp; e, ux) + sWx(x p ; e, v x ) + s 2 W 2 (xp; e, v t ) + • • • . (4.2) 

Substituting Eq. H4.2[) into the boundary conditions (13.71) and (I3.8p . and noting 
that s = vx, we obtain 

2a lr*3 U l» + a lr3a W 0,a + M^ a lra3 U 2a + a\ T 3a W l,cc) 

+ -ji^lr^U^Uxa + 4:a 2 3ra333 W U la 
+Hr3a^U ^W , a + 2al r333a W W , a ) 

+ "g"( 6a 3ra3»£AK^la^e,7^0K,A + ^a 3Ta333XK U la W U 0K:X + 603^33333 U la W Q 
+ ^ a 3r3a^XK W 0,aUo^rjUo Ki X + 6a| r3a ^3 3 Wq^Uq^Wq + 80^3333 Wq^Wq) 

+ O(e 3 ,v 1 e) = (r=l,2), (4.3) 

a 33f3aUoa,{3 + 03333^0 + V\ («33/3a^la,/3 + 3a 3 33 3 Wl) 

f 2 2 2 2 \ 

+ ( a 33a/3^^ l9 7^0/3,a + 20.33^33 W [/oft« + 0-333333 ^0 ) 

+ "^"( a 33a/3r;CA K ^0/3,a^05, ?? ^0« : ,A + ^33a/3rjC33 U 0P,» U 0^riW 
+3 a 33a/33333^0/3,aW / 2 + a 33333333^0 ) 

+ O(e 3 ,v 1 e) = 0. (4.4) 

It should be noted that boundary conditions (14.31) and (14.4j) can be expanded 
to any needed order. But, for the problem we consider here, we omit 0(e 3 , v x e) 
and higher-order terms, since the purpose is to deduce the leading-order be- 
havior. We note that the above three equations contain eight unknowns Uo a , 
Wo, U Xa , W x and U 2a . Thus we need another five equations to have a closed 
system. 

Substituting Eq. (14. 2p into (13. 5p . the left-hand side becomes a series in s, 
and all the coefficients of s n (n = 0, 1, 2, • ■ ■) should be zero. As a result, we 
have two sets (r = 1, 2) of infinitely many equations. Among them, we only 
consider those contain the eight unknowns as in (14.31) and (14.41) . Actually, from 
the coefficients of s° and s 1 , we obtain 

a L 7 /3^o/3,Q 7 + 2a 3T 3 a L r i Q + (al ra3 + 0^33)^0,0 
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+ e ( a aT 1 f3 V iUoi3, ai Uo^ v + 2a 3ra3vi Uo^ ri U la + a aT1 p 33 U p tCr/ W 

+ ( a 3r3/3r?5 + a |r33^) U^Wq^ + 2al ra333 Ui a W 
+ ( a 3r3/333 + a /3r3333) W W 0t p) 

+ e 2 (H 5 ) = (r = l,2), (4.5) 



a otryP U W,on + 12a 3r3a^2 a + 3( a 3ra3 + a ar33)Wl," 
+ e ( a /3r7a»?5(^ r i«,/97^0e,7 ? + U OafaU + ^(a^g^g + a^p^Ut^pUxr, 
+ l 2a 3rQ3» ? ^0C,r ) ^ r 2a + °^ T 7a33^'la,/97 1 ^0 + (3 a 3r3^ + a /3r33r^) ^fl^O,/? 
+ 12aL3333^2aWb + 3(^33 + ^3333) (W W ltf} + WjWqJ) 
+ 2 a 3Ta3r 1 tiUu,r 1 Ula + 2a^ T ^ 3 L r i ?? W / o,a / 3 + a ar3/ 33^Wo,a/?Wo,»7 

+3a^33^,« J ^i + ISa^^iaWi + (9aj r3a33 + 3^3333)^0,^1 

+ e 2 (H 6 ) = (r=l,2). (4.6) 

Similarly, substituting Eq. (14.21) into (I3.6p . we have a set of infinitely-many 
equations. We only use the equation coming from the coefficient of s° since 
only it contains the eight unknowns mentioned before. The equation takes the 
form: 



2 ( a 33,3a + a a3/33)^l«,/3 + a a3/33^0,Q/3 + 603333^1 

{ i 2 2 \ / 2 2 \ 

+ e ( 2 ( a a3/337?€ + a 33/3ar,5)^la,/3^,r; + 2{d a3 £ 333 + a.33^33) U la £ W 

+2a^3a ? ?C3 f/ oa,r ? /3^ + 4a| 33a g3 W , a ^i? + 4a^ 3 33 ?3 L r i Q L r i § 

+ a a3/33^ W/ 0,a/3^ / "0C,r, + ^37035 ^0a,/9 7 ^0,£ + ^^333^0,0/9^0 

+ 2aL 3 33e(t / la + WV)WV + 6a^33^C/ ^ a Wi + 6CI333333WW1) 

+ e 2 (H 7 ) = 0. (4.7) 

Now, the equations H4.3[) to (14. 7p provide the eight governing equations for 
the eight unknowns Uo a , Wo, Ui a , W\ and JJ-ia and we have a closed system 
to work with. 



To further simplify the two-dimensional system of Eqs. (14. 3[) to (14.71) . we shall 
further use the smallness of the parameter e through asymptotic expansions. 
By a perturbation method, from fl4.5[) we obtain 
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+ (^3(«L33/3A + a ta63f3x) ~ ( a 3 7 53 + ^S^^a^xWoX^Wo^ 
+ (A3a>0 a Sy£r) ~ a /3C7<5 a 3a3C€r))^05,/37^0r)^] 

+ e 2 (H 8 ), (4.8) 

where Ai and B{ (i — 1, 2, • • ■) are defined in Appendix A. Substituting (14.81) 
into (14.71) . we obtain 



Wl = it 3 W °' aa + ^X^ a W«Uo a ,w + 0(e), (4.9) 

where we only give the expression of the leading-order term of W\ since the 
higher-order terms have no influence on the final asymptotic equations. Sub- 
stituting (Q|) and (l4~9l) into fl4~6l) . we obtain 



-A X (A 2 + A 3 )(A 2 W ,a^ - aJatfHW)) + O(e). (4.10) 
Substituting gSD to (OCT) into g3) and (Ojl . we obtain 



Ag^Or.aa + A 2 W , T + (A 2 + A 3 )U 0a , ar 
- j{A 3 U 0Tiaa ^ + {3A 2 + 2A 3 )W , a , 

+ ^((A 3 a% 1&33 - B 4 al TlS )U 0S ,frW + (Ai(B 7 - B 4 ) + A 3 B 2 )W W , 7 



+ {3A 2 + 2A 3 )W ,aar + 3(A 2 + A 3 )U Qa ^ T ) 

e 

A 3 

/ 2 1 /i2 1 2 \ 

+ ( a 3Ta3^ a /3aA 7 + A 3 <Xp TXl ^ — 0, p^a^^^Uox^Uo^^ 

.2 „2 \ A „2 



+ (^l( a 3ra3/3A _ a 3r3a(3x) ~ ^3^33^) Wo,a) 

+ e 2 (tf 9 ) = 0, (4.11) 

^Oa.a + ^Wq + y (AtW/3/3 + A 2 W , Qa ) 
+ 9( a 33a/3^^05,»?^0/3,a + 2B 2 W Uo at a + B\Wq) 



2 
f _ 2 

6 

2 i r< tj/3\ 



+3C 2 C/ a,aW Q 2 + dWtf) = 0, (4.12) 

where we have omitted the terms higher than 0(e 2 , v{). 

For the two side surfaces, we suppose that they are traction-free. Also, since 
for a thin plate they are much smaller than the top and bottom surfaces, for 
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the latter we need the boundary conditions to be satisfied at every point (cf. 
( 14. 3 P and (I4.4[) ) while for the former we only require the boundary conditions 
to be satisfied in an average sense along the thickness. By integrating the 
traction-free boundary conditions E 2 j = at x 2 = (where v 2 = b 2 /l 2 ) 

along the thickness from to a, we obtain 

- ^{2(A 2 + A 3 ){W + U iA + 1/02,2X12 + A 3 (U Q1 , 2 + U 02jl ) iaa ) 
o 

+ 2( a 21a/37<5^0/3,Q^05, 7 + ^ a ll33a/3Uop,aW ) 

+ ~^( a 21al3-rSK\Uof3,aUoS n Uo\ jK + 3a 2 i33 aj g 7( j?/oa,a^05,7 Wo 

+ 3^ 13 33 3a/ 3?7oa I aWo 2 ))U a =±V^ = 0, (4.13) 



2^02,2 + ^2(^0 + ^) 

- ^(2(A 2 + A 3 ){W + U X,1 + U 02 , 2 ),22 + (^^02,2 + A 2 (l^ + C/ l,l)),aa) 


+ ^(^a^^OAa^^ + 2a 22a/333 ^0/3,0 Wo + -^Wq 2 ) 
e 2 

+3^3333^^ + 3CW 3 )U =±v ^ = 0. (4.14) 



^2^0,2 + (A 2 + A 3 )[/ a,a2 + A 3 U 02 , aa 
- ^((3A 2 + 2A 3 )W ,aa2 + 3(A 2 + A 3 )f/o a ,a/3/32 + A 3 U 02 , aa ^) 

+ e(# 10 ) + e 2 (#ii)|, 2=±v ^ = 0. (4.15) 

Eqs. (14. lip and (14.121) are the three asymptotically-valid governing equations 
for the three unknowns U 0a and W , among which U 0a are the two displace- 
ment components of a point in the middle plane and Wq is the normal strain 
(along the thickness direction) of that point. Since the two-dimensional system 
of Eqs. ( 14.111) and ( I4.12p together with the six boundary conditions ( 14. 131) to 
( I4.15P are derived from the three-dimensional field equations, once the solu- 
tion of this system is obtained, the three-dimensional displacement field (thus, 
also the strain and stress fields) can be easily calculated. 



Remark: It is easy to see that these governing equations and boundary condi- 
tions have two significant features which are different from the equations for a 
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standard plane-stress problem. Firstly, the out-plane normal strain Wo is cou- 
pled with the in-plane displacement components U 0a . And secondly, there are 
some v\ terms and the orders of the derivatives of these terms are two-order 
higher than the other terms, which indicate the influence of the thickness of 
the plate. We shall see later that the thickness has an important influence 
on the bifurcations. Thus, a model based on a plane-stress problem may be 
defective for capturing the instability phenomena in a thin plate. 



5 Asymptotic Normal Form Equations for a Thin Strip 

Now we consider the case that the plate is a thin strip in terms that both the 
thickness and the width of the plate are much smaller than the length (this is 
in agreement with of the experimental setting of Shaw and Kyriakides 1998). 
Thus, besides v\ being small, v 2 is also small. As a result, —^/T>2 < x 2 < ^/v^ 
is a small variable. From Eqs. fl4.lip - fl4.15l) . it is clear that the unknowns are 
functions of the variable X\[=: x), the small variable x 2 (=: y) and the three 
small parameters e, z^i, i/ 2 , i.e., 

Uoa = U a{x,y,e, vi,is 2 ), W = W (x,y;e,u 1 ,u 2 ). (5.1) 

We assume that the unknowns are sufficiently smooth in y and seek the 
series expansions in the small variable y: 

U 01 =u (x) + y 2 u 2 {x) + y 4 u 4 (x) + y & u % {x) H , 

+v^2 y ■ {u\{x) + y 2 u 3 (x) + y 4 u 5 (x) + y 6 u 7 (x) H ), 

U 02 = Jv~ 2 ■ (v (x) + y 2 v 2 (x) + y A v A (x) + y 6 v 6 (x) H ) 

+y ■ (vi(x) + y 2 v 3 (x) + y\(x) + y 6 v 7 (x) H ), 

W = wo(x) + y 2 w 2 {x) + y 4 w 4 (x) + y 6 w 6 (x) H 

+v / ^2 V ■ {wi{x) + y 2 w 3 (x) + y 4 w 5 (x) + y 6 w 7 (x) H ), (5.2) 

where ^/T^ is introduced into the expansions based on the assumption that the 
maximum non-dimensional lateral displacement of a point in the center line of 
the middle plane is 0( v /z^) (i.e., Uo 2 \ y =o = 0{ y /u 2 )) 1 since in the experiment 
there was only a small bending (see Shaw and Kyriakides 1998). 

Substituting f!5.2p into (14. 13|) and (14.14D and omitting terms higher than 
0(e 2 , v 2 ) and with some manipulations, we obtain 

M{u\ + vox + ^ 2 (3m 3 + v 2x )) 
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- — (2A 2 (2v 2x + w lx + u lxx ) + A 3 (6u 3 + 6v 2x + 2w lx + 3u lxx + v 0xxx )) 
o 

+ t(ui(B A (u 0x + vi) + B 5 w ) + v 0x (B 7 (u 0x + vi) + B 8 w )) 

+ -q(3(v 0x (C 12 uI x + C 12 vf + 2C 15 v 1 w + C 13 Wq + 2u 0x (C u v 1 + Ci 5 w )) 

+ui(C7 5 i4 e + C 5 v\ + 2C 8 t;iWo + Oo + 2u 0x (C 7Vl + C 8 u>o)))) = 0, 

(5.3) 

A 3 (2u 2 + v lx + z/ 2 (4w 4 + v 3x )) 

- ^-(AA 2 (3v 3x + w 2x + u 2xx ) + A 3 (24w 4 + \%v 3x + 4w 2x + Qu 2xx + ^i^)) 
o 

+ e(2u 2 (B i (u Qx + ui) + .65^0) + v lx (B 7 (u 0x + f 1) + B 8 w )) 

+ 7r(3(wix(C , i2Mox + Ci2^i + + C 13 Wq + 2u Qx (C u v 1 + Ci 5 u> )) 

o 

+2u 2 {C 5 u 2 0x + C b v\ + 2C 8 v lWo + C 6 w 2 + 2u 0x {C 7Vl + C 8 w )))) = 0, 

(5.4) 

A 2 {u 0x + w + v l ) + 2A 3Vl + v 2 {A 2 [u 2x + w 2 ) + 3{A 2 + 2A 3 )v 3 ) 

- ^(6(3A 2 + 4^3)^3 + 2(3A 2 + 2A 3 )(w 2 + u 2x ) 


+A 1 v lxx + A 2 (u;oxx + u Qxxx )) 
+ e -(B 2 u\ x + B x v\ + 25 2 ^iw + B 2 iug + 2uv x (B 2 v x + 5 3 w )) 

e 2 

+ — {C 2 u\ x + + 3C 2 v\w + 3C 3 v x wl + C 2 wl 
o 

+3ug a .(C 3 Ui + C4W0) + 3u 0x {C 2 v 2 + 2C^V\Wq + C A wl)) = 0, (5.5) 

^2 (u lx + w x + 2v 2 ) + 4A 3 w 2 + is 2 (A 2 (u 3x + Av 4 + w 3 ) + 8A 3 v A ) 

- ^(24(3A 2 + AA 3 )v 4 + 6(3A 2 + 2A 3 )(u> 3 + u 3x ) 


+2A 1 v 2xx + A 2 (w lxx + u lxxx )) 
+ e(2u 2 (B 4 ui + B 7 v 0x ) + (B 7 ui + B 4 v 0x )v lx + 2B 2 v 2 (u 0x + w ) 
+u lx (B 2 u 0x + B 3 w ) + {B 3 u 0x + B 2 w )wi 
+v 1 (2B 1 v 2 + B 2 (u lx + w 1 ))) 

e 2 

+ — (3{C 2 ul x ui x + ^Ci 4 u 0x u 2 vq x + 2C 3 u\ x v 2 + 2C 4 uo x u lx wo 
o 

+4:Ci 5 u 2 v 0x w + ACiUQ X v 2 WQ + C^uixWq + 2C 3 v 2 u>o 
+2i>o x t>ix(C7Uox + C 8 w ) + 2u 1 (vi x (Cuu 0x + C 12 vi + C 15 w ) 
+2u 2 (C 7 u 0x + C 5 v t + C 8 w )) + C^u^Wi 
+2C4M0XW0W1 + C 2 wlwi + 2^i (^(2(712^2 + Cs^ix) 
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+2C 2 v 2 (u 0x + w ) + u lx (C 3 u 0x + C 4 w ) + {C 4 u 0x + C 3 w )wi) 
+v 2 l (2C l v 2 + C 2 (u lx + w l )))) = 0. (5.6) 



We note that the above four equations contain fourteen unknowns u — u 4 , 
Vq — v 4 , w — w 3 and to have a closed system, we need another ten equations. 

Substituting (15.21) into Eq. (14.111) for r = 1, all the coefficients of y n 
(n = 1, 2, • ■ ■) should be zero. So, we have a set of infinitely-many equations. 
However, only the coefficients of y°, y 1 and y 2 contain the fourteen unknowns 
mentioned above, which yield the following three equations: 

A 2 (u 0xx + w 0x + v lx ) + A 3 (2u 0xx + v lx + 2u 2 ) 

„ \- rl 2\' JU '[) xxxx + §u 2xx + 3wq xxx + Qw 2x + 18t>3 X . + 3v\ xxx ) 
o 

+A 3 (4u 0xxxx + 10u 2xx + 2Au 4 + 2w 0xxx + Aw 2x + l%v 3x + 3v lxxx )) 
+e(B 2 u 0x v lx + B 7 u 0x v lx + B 3 v lx w + B 8 v lx w + Uq xx (BiUq x + B 2 w ) 
+2u 2 (B 4 (u 0x + v{) + B 5 w ) + B 2 u 0x w 0x + B 2 w w 0x 
+v 1 (B 2 u 0xx + (B 2 + Bt)v Xx + B 3 w 0x )) 

+^(C 12 u 2 0x v lx + C 2 ul x v lx + 2C 15 u 0x v lx w 
+2C 4 u 0x vi x w + C 13 v lx wl + C A v lx wl 

+u 0xx (C 1 ul x + 2C 2 uq x wq + C 3 Wq) + 2u 2 (C 5 ul x + C 5 vf + 2C s viw 
+C e w 2 + 2u 0x (C 7 v 1 + C 8 w )) + C 2 ul x w 0x + 2C 3 u 0x w w 0x 
+C 2 w 2 .w 0x + v\(C 3 uq xx + (C12 + C 2 )v lx + C 4 w 0x ) + 2v 1 (u 0x (C 2 u 0xx 
+ (Ci4 + C 3 )v lx + C 4 w 0x ) + w ((C 15 + C 4 )v lx + C 4 (u 0xx + w 0x )))) = 0, 

(5.7) 

A 2 {u lxx + w lx + 2v 2x ) + 2A 3 (u lxx + v 2x + 3u 3 ) + 0(e) = 0, (5.8) 

A 2 (u 2xx + 3v 3x + w 2x ) + A 3 (2u 2xx + 3v 3x + 12u 4 ) + 0(e) = 0. (5.9) 

The 0(e) terms in (I5.8P and (I5.9P will not make contributions to the final 
results and are thus not written out. Similarly, substituting (15.21) into Eq. 
(14. lip for r = 2 and from the coefficients of y°, y 1 and y 2 , we obtain 

A 2 (u lx + u>i + 2v 2 ) + A 3 (u lx + v 0xx + Av 2 ) 

—^{A 2 {3u lxxx + l%u 3x + 3wi xx + 18w 3 + 72v 4 + 6v 2xx ) 


+A 3 (3u lxxx + 18u 3x + 2w lxx + 12w 3 + 96v 4 + I0v 2xx + v 0xxxx )) 
+e(2B 2 u 0x v 2 + 2B x viv 2 + 2B 2 v 2 w + v 0xx (B 4 (u 0x + vi) + B 5 w ) 
+u lx {{B 2 + B 7 )(u 0x + Vl ) + (B 3 + B 8 )w ) 
+v 0x {B 4 u 0xx + 2B 7 u 2 + 2B 4 vi x + B 5 w 0x ) + B 3 u 0x wi 
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+Ui(B 7 u 0xx + 2B 4 u 2 + 2B 7 v lx + B 8 w 0x ) + B 2 v 1 w 1 + B 2 wqWi) 
+—{C\ 2 ul x ux x + C 2 ul x u lx + 2C 5 uo x u 0xx v 0x + AC ia uq x u 2 vq x 

+4C 7 u 0x v 0x v lx + 2C 3 u 2 0x v 2 + 2C 15 u 0x u lx w + 2C 4 u 0x u lx wo 
+2C 8 u 0xx v 0x w + AC 15 u 2 v 0x w + 4C 8 v 0x v lx w + 4:C 4 u 0x v 2 w 
+C 13 u 1x Wq + C A ui x wl + 2C 3 v 2 Wq + C 4 ul x wi 
+v 0xx (C 5 ul x + 2C 8 u 0x w + C & wl) + 2C 8 u 0x v 0x w 0x + 2C e v 0x w w 0x 
+2u 1 (2((C u u 0x + C 12 vi)v lx + C 15 v lx w ) + u 0xx (C 12 u 0x + C u vi + C 15 w ) 
+2u 2 (C 7 u 0x + C 5 vt + C 8 w ) + C 15 u 0x w 0x + CmViWox + C 13 w w 0x ) 
+2C i u 0x w w 1 + C 2 wlw 1 + v\{{C 12 + C 2 )u lx + C 5 v 0xx + 2C x v 2 + C 2 w x ) 
+2v x (2C 2 v 2 {uq x + w ) + v 0x (C 7 u Qxx + 2C 12 u 2 + 2C 5 Vi x + C 8 w 0x ) 

+^o((Cl5 + Ci)ui x + C 8 V 0xx + C3W1) 

+u 0x ((C u + C 3 )u lx + C 7 v 0xx + C 4 wi))) = 0, (5.10) 
2A 2 [u 2x + w 2 + 3v 3 ) + A 3 (2u 2x + v lxx + 12v 3 ) + 0(e) = 0, (5.11) 

3A 2 (u 3x + w 3 + Av A ) + A 3 (3u 3x + v 2xx + 24v 4 ) + 0(e) = 0. (5.12) 

Substituting (15.21) into (14.121) and from the coefficients of y°, y l , y 2 and y 3 , we 
obtain 

A 2 (u 0x + w + v 1 ) + 2A 3 w 

—pr{A 2 {3u Qxxx + 6u 2x + 3w 0xx + 6w 2 + 18v 3 + 3v lxx ) 
6 

+A 3 (6u 0xxx + I2u 2x + I8v 3 + 6v lxx )) 

+ ^{B 2 u 2 0x + B 2 v\ + 2B 2 v lWo + B lW 2 + 2u 0x {B 3Vl + B 2 w )) 

+—(C 2 ul x + C 2 v\ + 3C 3 v 2 w + 3C 2 wi^o + c i w o 
6 

+3u 2 0x {C i v 1 + C 3 w ) + 3uo*(C 4 v? + 2C A V!W + C 2 v%)) = 0, (5.13) 
A 2 {u\ x + wi + 2v 2 ) + 2v4 3 wi 

— jr(A 2 (3ui xxx + 18u 3x + 3w lxx + 18w 3 + 72w 4 + Qv 2xx ) 
+A 3 {6u lxxx + 36u 3x + UAv 4 + I2v 2xx )) 

+e(2u 2 (B 5 u x + B 8 v 0x ) + {B 8 ui + B 5 v 0x )v lx + B 2 u lx (u 0x + w ) 
+2v 2 (B 3 u 0x + .BaWo) + (B 2 u 0x + B 1 w )w 1 + Vi(B 3 u lx + 2B 2 v 2 + B 2 w{)) 

+ — {3{C 2 ul x u lx + ^Ci 5 u 0x u 2 v 0x + 2C 3 u 0x Ui x w + 4C 13 u 2 v 0x w 
o 

+2C^u 2 0x v 2 + 4:C 4 u 0x v 2 Wo + C 2 (u lx + 2v 2 )w 2 Q + 2v 0x v lx (C 8 u 0x + C 6 w ) 
+2ui(v lx (C 15 (u 0x + vx) + C 13 w ) + 2u 2 (C 8 (u 0x + vx) + C 6 w )) 
+C 3 ul x wi + 2C 2 u 0x w wi + Ciwlwi + v\{C^u\ x + 2C 2 f2 + C 3 w\) 
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+2v 1 (v 0x (2C 15 u 2 + C s v lx ) + C 4 u lx (u 0x + w ) + 2v 2 (C 4 u 0x + C 3 w ) 

+ (C 4 u 0x + C 2 w )w 1 ))) = 0, (5.14) 

A 2 (u 2x + w 2 + 3v 3 ) + 2A 3 w 2 + 0(e) = 0, (5.15) 

A 2 (u 3x + w 3 + Av 4 ) + 2A 3 w 3 + 0(e) = 0. (5.16) 



Now, we have a closed system: the fourteen equations (15.31) to (15.161) provide 
the governing equations for fourteen unknowns u - w 4 , v - t> 4 and w - w 3 . 
Please note that the boundary conditions (I4.15P can be satisfied automatically, 
and so it will not be used. It is still formidable to analyze the solutions of this 
system directly. Next, we shall use the smallness of the three parameters e, V\ 
and v 2 to proceed further. 

By a regular perturbation expansion, from (I5.3P we can obtain the expression 
of u\ as a function of uq and other unknowns. By substituting this expression 
of Mi into the remain equations of (15. 3p to (15.161) . we can eliminate u\ from 
all these equations. Similarly, by solving the resulting equation of (15.71) xAo 
+ (15.41) x (7A> + 4:A 3 ), we can obtain u 2 as a function of u , v , . . .. And by 
substituting this expression of u 2 into the remain equations of (15. 3p to (I5.16p . 
we can eliminate u 2 from all these equations. Similarly, we can express Ui-u 4 , 
Wq-w 3 and V\-v 4 in terms of Uq and vq. The concrete forms are 



A 2 v x - 3{3A 2 + 2A 3 )v 2 

ui — —v 0x -\ —— — ■ — — v 0xxx 

6{A 2 + A 3 ) 

+ e 2 (A 2 + + A 3 ) U0xV ° x + e2 ^ Hl2 ^ v °^ ( 5 - 17 ) 
Aj 

U2 ~8(A 2 + A 3 )(2A 2 + A 3 ) U0xx 



A 2 {3A 2 + 2A 3 ) 



192(A 2 + A 3 ) 2 (2A 2 + A 



3/ 



■((15 A 2 + 8A 3 )^ - 3(3A 2 + 2A 3 )u 2 )u 0x 



e(3A2 + 2A3) -{AAlB 2 -AAl{A 3 + B,) 



32A 3 (A 2 + A 3 f(2A 2 + A 3 
-2A 2 A\{AB % -B 2 + B 3 ) - A*A 3 {bA 3 + 4(5 6 - B 2 + B 3 )))u 0x u 0xx 
+ e 2 (H 13 )ul x u 0xx , (5.18) 
5-42 + 4A 3 

U3= l2(A 2 + A 3 ) V0xxx + O(e) > (5 ' 19) 



A 2 (19A 2 + 20A 2 A 3 + 4A 2 ^ 
384(A 2 + A 3 ) 2 (2A 2 + A 3 ) 



u 4 = — OOA , 4 — , , x 9 / , — — rr«0xxH + 0(6), (5.20) 
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32{A 2 + 2A,)\2A 2 + A 3 )™ ' (5 ' 22) 



A 2 A 2 (5(3A 2 + 2A 3 )v 2 - A 2 v l ) 

Wl ~ 2{A 2 + A 3 f 0xx + 24(A 2 + A 3 ) 2 

A 2 (19A 2 2 + 26A 2 A 3 + 82A 2 3 ) 
" 6 lQ{A 2 + A 3 f{2A 2 + A 3 ) UoxxVox 

AA 2 B 2 + A\{B X + 3B 2 - 4E 3 ) + 4A 2 A 3 (£ 2 - B 3 ) 

KM + a 3 ) 3 

A 2 (3A 2 + 2A 3 ) 



+ ei 



-)«0a^0 a 



4(A 2 + A3) 2 

+ e 2 (ifi 4 U0xxl'0s + #15^0^0x0:) ^(te, (5.23) 

3A 2 (3A 2 + 2A 3 )z/ 2 - (71A 3 2 + 208A 2 A 3 + 180A 2 A 2 + A8A\)v 1 

+ 96(A 2 + A 3 ) 2 (2A 2 + A 3 ) U ° xxx 
_ 4A\B 2 + A 2 (Ex + 3£? 2 - 4£? 3 ) + 4A 2 A 3 (B 2 - B 3 ) 2 

6 16(A 2 + A 3 ) 3 Mte 

+ 6 2 (iJ 16 )^, (5.24) 



^ = - 24 2 (1 2 2 + aV O ^ + 0(£) ' (5 ' 25) 
A 2 (A 2 2 + 8A 2 A 3 + 4Ap 

96(A 2 + A 3 ) 2 (2A 2 + A 3 ) 

A 2 (3A 2 + 2A 3 )((7A 2 + 2^3>^Vi) 

U2 4(^2 + ^)^ 48(A 2 + A 3 ) 2 V ° xxxx 

A 2 (3A 2 + 2A 3 ) 2 



+ e 



32(A 2 + A 3 ) 2 (2A 2 + A 3 ) 
AA\B 2 + A|(.Si + 3ff 2 - 4B 3 ) + 4A 2 A 3 (£ 2 - B 3 
[ 16(A 2 + A 3 y 

A 2 (3A 2 + 2A 3 ) 



-)u xVo x 



8(A 2 + A 3 ) 2 

+ e 2 (H 17 u 0xx v 0x + Hi 8 u 0x v 0xx )u 0x , (5.27) 

Vi = -2(A 2 A Ia 3 ) U0x 

A 2 ((71A 2 + 98A 2 A 3 + 32A 2 )u 1 - 3(A 2 + 2A 3 )(3A 2 + 2A 3 )u 2 ) 
+ 96(A 2 + A 3 ) 2 (2A 2 + A 3 ) 
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AA 2 B 2 + A 2 {B X + 3B 2 - AB 3 ) + AA 2 A 3 {B 2 - B 



— e 7—. r-^ u 



O.r 



16(A 2 + A 3) 

+ e 2 (H 19 )u 3 0x . (5.28) 

By inserting (15. 1 T[) to (15.281) into (15.71) and (I5.10p . we obtain two governing 
equations for the two basic unknowns Uq and Vq\ 



uqxx + 2D 1 eu 0x uo xx + 3D 2 e 2 ul x u 0xx - ^u 0xxxx = 0, (5.29) 



Ev 

E((eu 0x + {Di - l)e 2 ul x )v 0x ) x —v 0xxxx = 0, (5.30) 



where 



A 3 (3A 2 + 2A 3 ) 

E= \ \ — - 5.31 

A 2 + A 3 K J 

is the Young's modulus, and 



Di = — r— ; ; ; r^UAZBx + 12A 2 A 2 JB 1 - B 2 ) 

8A 3 (A 2 + A 3 ) 2 (3A 2 + 2A 3 y 3 2 31 1 v 

+ QA 2 2 A 3 (2B 1 - 3B 2 + B 3 ) + 3A\(B 1 - 3B 2 + 2B 3 )), (5.32) 

D 2 = I 

96A 3 (A 2 + A 3 )\3A 2 + 2A 3 ) 

■ {{1QA\{A 3 C X - 3B 2 ) + 16A 2 Al(-6B 2 2 + 6B 2 B 3 + 5A 3 d - AA 3 C 2 ) 
+ 18^(d - 4C 2 + 3C 3 ) + 8A 2 2 A 2 3 (-3B 1 B 2 - \hB 2 2 + 24B 2 B 3 - 6B 2 
+ 20A 3 d - 32A 3 C 2 + QA 3 C 3 + QA 3 C 4 ) + 8A 3 A 3 (-35 1 fi 2 - 9B 2 
+ 3B Y B 3 + 21B 2 B 3 - \2B\ + 204A - 50A 3 C 2 + 18A 3 C 3 + 12A 3 C 4 ) 
+ A\(-3Bl - 18£i£ 2 - 27B 2 + 24^1^3 + 72^2^3 - 485| + 82A 3 d 
- 280A 3 C 2 + 150A 3 C 3 + 48A 3 C 4 )), (5.33) 
2(3A 2 + 2A 3 )u 1 

v= K . : - A { . 5.34 
3(A 2 + A 3 ) 1 ; 

It can be seen that Di and D 2 are constants which depend on material con- 
stants and v is proportional to the small parameter V\. As U\ — u 4 , V\ — v 4 , 
w — w 3 are expressed in terms of u 0x and v 0x , once u 0x and v 0x are found and 
all these quantities can also be found. 

Integrating (I5.29j) with respect to x once, we obtain 



u 0x + D x eu\ x + D 2 e 2 u\ x - -« « = C, (5.35) 



19 



where C is an integration constant. It is important to find the physical meaning 
of C, since we aim to investigate the instability phenomena as the physical 
parameters vary For that purpose, we consider the resultant axial force T 
acting on the material cross section that is planar and perpendicular to the 
strip axis in the reference configuration, and the formula is 



b a 

T = J J E n dX 3 dX 2 . (5.36) 

— b —a 

By using Eqs. (I2^jl . (4.2), (4.8)-(4.10) and (PTjl - rtQg]) in (15T36D . it is possible 
to express En in terms of Uq x and Vq x - Then, carrying out the integration in 
f )5.36p . we find that 



T = 4abEe{u 0x + D x eul x + Aj^ujL - ^u 0xxx ). (5.37) 

Comparing Eqs. (15.351) and (I5.37p . we have C = .T E . Thus, we can rewrite 
( ESSJ) as 



eu 0x + Di(euQ X ) 2 + D 2 (eu 0x ) 3 - ^-eu 0xxx = -r^-=. (5.38) 

If we retain the original dimensional variable and let U = uqx = £Uo x (where 
X = Ix = X\), we have 



a 2 

U + D 1 U 2 + D 2 U 3 -^U X x=l, (5.39) 



where 



^^E (5 - 4 °» 
is the engineering stress and a\ = ^fp~^\^ • 



Similarly, integrating ( 15 .301) with respect to x once, we obtain 



E[eu 0x + (£>i - ^e 2 ^^,]^ —v 0xxx = D, (5.41) 

where D is an integration constant. To find out the physical meaning of D, we 
consider the resultant shear force Q acting on the material cross section that 
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is planar and perpendicular to the strip axis in the reference configuration, 
and the formula is 



b a 



Q = J J E 12 dX 3 dX 2 . (5.42) 



—6 —a 



By using Eqs. (l2T4"j) . (4.2), (4.8)-(4.10) and (l5TT71) - (l5^8|) in §M>, it is possible 
to express £12 in terms of Uq x and Vq x . Then, carrying out the integration in 
(ET35]) . we find that 



Q = 4:abE,fv2~e((eu 0x + (D 1 - l)e 2 ul x )v 0x - —v 0xxx ). (5.43) 

Comparing Eqs. (15.411) and (15.431) . we have D = 4ab < ^ € ■ Thus, we can rewrite 
(ESP as 



AabE((eu 0x + (D l - l)e 2 u 2 0x )ey/i^v 0x - —e^v 0xxx ) = Q. (5.44) 



If we retain the original dimensional variables and let V = ey^^oa 
(= M 2,i|(x 2 ,x 3 )=(o,o), where m 2j i is defined in (2.2)), we have 



EA(U+(D 1 -1)U 2 )V-EJV XX = Q, (5.45) 

where A is the area of the cross section and J is the moment inertia around 
the X3-axis of the cross section. 

Once U (i.e., u 0x ) and V (i.e., v 0x ) are found from equations f!5.39j) and 
(I5.45p . all the other physical quantities can be calculated immediately. Also, 
since these two equations are derived in a mathematically consistent manner 
and contain all the required terms to yield the leading-term behavior of the 
original three-dimensional problem, we call them to be the asymptotic normal 
form equations of the governing three-dimensional nonlinear PDE's (2.5) with 
the traction-free boundary conditions on the top/bottom and two side surfaces 
under the given end axial resultant T and end shear resultant Q. 



Remarks: The results established by us can also provide some useful infor- 
mation on some classical results obtained before in an ad hoc manner. The 
details are discussed below. 

In the case that U is constant, and the nonlinear terms in (15.391) are ne- 
glected, we have U = 7. If the shear force Q is zero and the nonlinear term 
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U 2 is neglected, then Eq. (15.451) equation becomes the linearized Euler buck- 
ing equation, which was obtained based on the hypothesis that the bending 
moment is proportional to the curvature. Here, we have derived this classical 
equation based on a three-dimensional setting in a mathematically consistent 
manner without such a hypothesis. 

We can also obtain the moment M acting on the X 3 -axis of the material 
cross section. The formula is 



b a 

M 



- J J X n (X 2 + u 2 )dX 3 dX 2 . (5.46) 



-b -a 



where u 2 is the displacement component along the X2-axis direction. By the 
same procedure as above, carrying out the integration in (I5.46p . we find that 



M = -4ablEe^2-((eu 0x + iVOo - jV 0xx ). (5.47) 



If we retain the original dimensional variable and let V = It^fv^VQ 
(= u 2 \(x 2 ,x 3 )=(o,o)), we have 



M = EJVxx - EA(U + D X U 2 )V. (5.48) 
Clearly, we have 



V = V x - (5.49) 

Note that if we let M be zero and do not consider the effect of the nonlinear 
term U 2 , equation (15.481) also becomes the classical Euler bucking equation. 

It should also be noted that in general that Mx ^ —Q as can be seen 
from (I5.47P and (I5.43p . But in the case that Uq x is constant, denoting x the 
deformed coordinate of the (namely, x = X + uo), then by neglecting 

the (^(eV^ 2 , e 2 ^ 2 ) terms, we have 



f = = _Q. ,5.50) 

ax 1 + euox 

Thus, here we have deduced the restriction under which the above classical 
result for a beam is valid. 
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6 Solutions for an Infinitely-long Strip 



According to the experiments (Shaw and Kyriakides 1998, Sun et al 2000, 
Tse and Sun 2000), for the phase transitions in a thin strip due to ten- 
sion/extention, there are at least three instability phenomena: (i) there is a 
formation of the transformation fronts (manifested as a neck); (ii) the neigh- 
boring two transformation fronts incline a same angle with the axial axis; (iii) 
the transformation front can switch to an orientation with the opposite angle. 
We refer (i) as a necking-type instability, (ii) as a shear instability and (iii) as 
a front-orientation instability. 

Now, we analyze the asymptotic normal form equations (I5.39P and (15.451) 
in order to shed insight into the instability phenomena observed in the exper- 
iments. 

As mentioned in Section 2, we consider this class of non-convex strain energy 
functions such that in a one-dimensional stress setting with a homogeneous 
strain state the engineering stress-strain curve has a local maximum and a local 
minimum, which is same as that considered in the classical paper by Ericksen 
(1975). To the third-order material nonlinearity, the stress-strain relation in 
this setting is provided by 

U + DiU 2 + D 2 U 3 = -f 

(i.e., in (15.391) by setting Uxx = 0, since U is independent of A in a ho- 
mogeneous strain state). The requirement that the 7 — U curve has a local 
maximum and minimum is equivalent to 

D x < 0, D 2 > 0, 3D 2 < D\ < AD 2 . (6.1) 

The peak stress value 72, the valley stress value 71 and the Maxwell stress 
value 7 m can be expressed in terms of Di and D 2 (cf. Dai & Cai (2006)): 

2Df - 2{D\ - ?>D 2 fl 2 - 9DiD 2 

71 ~ rw\ ' 

2D\ + 2{D{ - 3D 2 ) 3 / 2 - 9L>iD 2 

72 ~ TW\ ' 

2D\ - 9D X D 2 

7m = 21D\ ■ (6 ' 2) 

The first normal form equation (I5.39P has the same form as that derived for a 
slender cylinder composed of an incompressible hyperelastic material (see Cai 
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and Dai 2006 and Dai and Cai 2006). In those two papers, it has been shown 
that this equation can be used to describe the necking-type instability and to 
capture the main features of the structure response (engineering stress-strain) 
curve. Thus, here we shall not study this equation and discuss the necking- 
type instability further. Instead, we concentrate on the other two instabilities: 
the shear instability and the front-orientation instability. For that purpose, we 
shall conduct a detailed study on the second normal form equation (15.451) . 

We focus on the case that there are two transformation fronts. If they are 
some distance away from the two ends of the strip and any other transforma- 
tion front (if present), without loss of generality, we can take the strip to be 
infinitely-long. Also we consider the case that the resultant shear force is zero 
(for an infinitely- long strip, it has to be zero otherwise the moment is infinite). 
Then to determine the solutions of the normal form equation (I5.45P becomes 
an eigenvalue problem with the eigenvalue equation 

EA(U + (£>! - l)U 2 )V - EJVxx = 0, 

or 

c 2 (-U + (1 - D^V + V xx = 0, (6.3) 



where c = ^JEA/EJ = ^ is a large parameter for a strip, and the boundary 
conditions are 



V = 0, at X = ±oo. (6.4) 

To solve (16.31) under (16.41) . one needs to solve the first normal form equation 
( I5.39P to get U(X;j). It can be viewed that 7 is the eigenvalue for (16. 3p . 
Denote /(X; 7) = -U + (1 - D^U 2 , then we can write fT3]) as 

V xx + c 2 f(X; 1 )V = 0. (6.5) 

Since c is a large parameter, it is possible to use the WKB method to construct 
the leading-order asymptotic solution (see Holmes 1998). There are three cases. 
If f(X; 7) > for —00 < X < +00, the general solution (to the leading order) 
is 



e ic J y/f(Xv)dX e -ic J y/f(X n )dX 

v=El +E2 imW' ■ (6 ' 6) 

If f(X; 7) < for —00 < X < +00, the general solution is 
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p c J y/-f{X;j)dX -c j y/-f(Xrr)dX 

V = E 3 r — r — A +E A - 



(6.7) 



If there are two turning points at Xq(> 0) and — X®, i.e., /(±Xq;7) = 0, 
f(X; 7) < for X < X < +oo and -oo < X < -X and f (X; 7) > for 
— X < X < Xq, the general solution is 



V = < 



for X > X + 61, 

C 2 Ai(C (X - X )) + C 3 Bi(C (X - X )) 

for Xe(X -5 1 ,X + 5 1 ], 
(/(X; 7 ))- 1/4 (C 4 sin(c/ X y/f(t- n )dt) + C 5 cos(c/ X y/f(P n )dt)) 

for Xef-Xo + ^Xo-^], 
C 6 Ai(C (X + X )) + C 7 Bi(C (X + X )) 

for X e[-X - 5 2 ,-X + 5 2 ), 
(-f(X; 7 ))- 1 /4(^ e - c l X x >/ = ?(^* + cr ge c /-x„ V 71 /^ 



for X < —X - S 2 , 
(6. 



where C±, C[, ■ ■ Cg are arbitrary constants, Ai(-) and Bi(-) are the Airy 
functions of the first and second kinds respectively, Co = \c 2 fx(X ; 7)| 1//3 and 
5i, 5 2 are quantities of 0(c~ 2 / 3 ). 

According to the experiments, the shear instability happens immediately 
after the formation of the two transformation fronts. Thus, we consider the 
anti-solitary wave solution of the first normal form equation (15.391) (which 
corresponding to the profile with two transformation fronts; cf. Figure 5 of 
Dai and Cai 2006). According to (6.9) of Dai and Cai (2006), the solution 
expression is 



U 



where 



92max ~ a 2 # 2 tanh (— ) 

1 -# 2 tanh 2 (^M 

z v 901 / 



(6.9) 



Ho 



92r, 



91 



a 2 - g x 



V2 



«2 -gi)(g 



2max 



9i)D 2 



(6.10) 
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and g\ is a double root of and g2max and a 2 are simple roots of 



l -U 2 + + \d 2 U a - jU - H = 0. (6.11) 

For this solution, it is easy to deduce that there are two turning points at ±X 
given by 



X = ga x arctanhW ^™ ax (6.12) 

V H 2 {a 2 - U ) 



where 



U = U(X on ) = T -^- (6.13) 

Thus, in this case the general solution is given by (16.81) . Upon using (16 .4p . it 
can be seen that C[ = C 8 = 0. Also, simple calculations show that fx(X ] 7 ) = 
U X {X Q]1 ), and thus C = \c 2 f x (X ; 7 )|V3 = \c?U x (X ; 7 )|V3. 



Remark: Since % is proportional to the strip thickness a (see the relation 
below (5.40)), thus the position of the turning point is also proportional to 
a. Thus, the thickness of the strip plays an important role in the instability 
phenomena. 

By using the matching conditions at the neighborhood of X = Xq, we can 
obtain relationships between constants C%, C 2 , C3, C4 and C5 in (16.81) . At the 
neighborhood of X = X , we have f(X; 7) ~ Ux(X ; j)(X — X ). The integral 
in (16. 8ft s can be written as 

X Xq Xq 

c J sjf(t;j)dt = c J y/f(t;i)dt-c J y/f(bi)dt 
x 

= /i(7)-^o(X -X)| 3 / 2 , (6.14) 

where 



Xo 



h(i) = c / v = a i c \ 



Qlmax 



2D? 



U{U - Uo 



dU 



at 



j^Pl k) + ~gU(Pl «) + e 2 K(k)) 



^ (a 2 - U)(g 2max - U) U - g x 

(6.15) 
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_ 9i(Uq-9i) a 2 (a 2 - Up) 

(« 2 - 9l)(g2max ~ 01 ) ' («2 - 52max)(«2 ~ 9l] 



VG(a>2 - g2max) _ _ ^^max ~ Up 

<?2^o( a 2 — Up)D 2 \ 92max{OL2 ~ Up 



^tanh 2 ^), ft = H 2 f3l (6.16) 
gai 

■) and /£"(-) are the complete elliptic integrals of the third and the first 
kind respectively. By comparing the asymptotic expansions of the second and 
the third forms in (16. 8p . we obtain 



C2 = ^1 Cj /4 |M*o; 7)r 1/4 P4(sin(A( 7 ) - cos(A( 7 ))), 

+C 5 (cos(/ 1 ( 7 )) + sin(/ 1 ( 7 )))], 
^3 = ^o 1/4 |^(X ; 7 )|- 1 / 4 [C 5 (cos(/ 1 ( 7 ))-sin(/ 1 ( 7 ))) ) 

+C 4 (cos(/ 1 ( 7 ))+sin(/ 1 ( 7 )))]. (6.17) 

By comparing the asymptotic expansions of the second and the first forms in 
we have 



^^-n-^C^lUxiXp^)^ 

C 3 = 0. (6.18) 

Similarly, by comparing the asymptotic expansions of the fourth and the third 
forms in (16. 8p . we obtain 



C G = yj\c^\U x {X Q ; 7 )|- 1 /4[_c 4 (sin(/ 1 ( 7 ) - cos(A( 7 ))), 
+C 5 (cos(A( 7 ))+sin(/ 1 ( 7 )))], 
|c 1/4 |f/x(X ; 7 )r 1/4 [C 5 (cos(/ 1 ( 7 ))-sin(/ 1 ( 7 ))), 
C 4 (cos(/ 1 ( 7 ))+sin(/ 1 ( 7 )))]. (6.19) 



By comparing the asymptotic expansions of the fifth and the fourth forms in 
( 16. 8p . we obtain 

c 8 = y^c 1/A \u x (Xp-n)\^c 6 , 

C 7 = 0. (6.20) 
Substituting (16.171) 9 into (16. 18|) q . we have 
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C 5 cos(/i( 7 ) + \) + C 4 cos(/ 1 ( 7 ) -j) = 0. (6.21) 
Substituting (16.191) 9 into (!6.2Up ^. we have 



C, cos(A( 7 ) + \) - C 4 cos(A( 7 ) - \) = 0. (6.22) 

For Eqs. (16.211) and (16.221) to have nontrivial solutions, there are two cases. 
One case is that 



C 4 = 0, /i( 7 ) = nvr + vr/4, n = 0, 1, 2, 3, • • • . (6.23) 

It is easy to see from (16.81) that in this case V is symmetric. Another case is 
that 



C 5 = 0, /x( 7 ) = mr - tt/4, n = 1,2,3, (6.24) 

and in this case K is anti-symmetric. 

Eqs. (I6.23P and (I6.24p are the eigenvalue equations for determining 7 . These 
two eigenvalue equations can be rewritten as a uniform expression 



(IN — lW , 
/i(7) = - ^ J -, N = 1,2,3,..-, (6.25) 

where odd iV represents the symmetric solution for V and the even N rep- 
resents the anti-symmetric solution. By substituting (16.151) into (16.251) . we 
obtain 



ei U((3l k) + «) + e 2 K{k) = - {2N 1)?r . (6.26) 

d\ 4 



It can be seen that the left hand side of equation (I6.26P implicitly depends on 
7 and the right hand side only depends on the width-thickness ratio and the 
wave number N. Thus, the width-thickness ratio is a key factor for determining 
the stress eigenvalues. 

Denote the eigenvalue corresponding to iV by 7e 7v- In the case that D\ = —18, 
D 2 = 100 and the Poisson's ratio (= 2{A f+ As) ) = 1/3, a = 0.01 and £ = 10 
(which means b ~ 0.1333), we obtain the first six eigenvalues as follows: 



7e i =lm + 6.045 x lO" 16 , 7e2 = 7m + 1.029 x 10" 41 , 
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7e3 = 1m + 1-751 x 10~ b( , 7e4 = lm + 2.980 X 10"^, 

7e5 = lm + 5.071 x 10- 119 , 7e6 = 7™ + 6.830 x 10~ 145 , (6.27) 

where the Maxwell stress 7 m = 0.0168. We note that these eigenvalus are very 
close. 

Denoting U and V the axial and lateral displacements of a point in the 
center line respectively, and without loss of generality letting the displacement 
at (0, 0, 0) be zero, we obtain 



U = J UdX, V = J VdX. 



(6.28) 



Curves of U, U, V, V for N — 1 and iV = 2 are plotted in Figure 1. It 
should be noted that the amplitude of the eigenfunction V (and then V) is 
undetermined. We have scaled all the amplitudes to unit for all variable in 
Figure 2. 
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Fig. 2. Curves of U, U, V, V as X varies. Left: N=l; Right: N=2. 
The coordinates in the current configuration are (up to the leading order) 



x = X + U - VY, 
y = Y + V- vUY, 
z = Z - vUZ. 



(6.29) 



The shapes of the thin strip corresponding to 7 e i to 7 e 6 are plotted in Figure 
3. The thickness of the strip in the current configuration is illustrated by flood 
contours. It should be noted that they represent six different modes and do 
not need to appear consecutively. 

From Figure 3, it can be seen that the transformation fronts are inclined 
with the strip axis. For N — 1, the inclined directions are the same and 
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Fig. 3. Shapes of the thin strip corresponding to j e \ to j e Q (from top to bottom). 
The thickness is illustrated by flood contours. 

the position of the turning point X is about 0.2675 (roughly speaking, X 
represents the position of the transformation front). For N = 2, the inclined 
directions are opposite and the position of the turning point X is about 
0.7619. In experiments it is found that, accompanying the formation of two 
phases, there is a stress drop (7 > j m ) and the phase fronts become inclined 
in the same direction (a shear instability). And then with the developing of 
the high-strain phase, the propagating phase front can become inclined in 
the opposite direction (orientation instability). Our analytical results seem to 
capture these features. 

The results obtained here shows that the inclination of the transformation 
front (shear instability) is a phenomenon of phase-transition-induced buckling. 
When the phase transition happens, there is a localized deformation (a neck). 
Due to that, the eigenvalue problem (6.5) has a turning point at Xq, which in 
turn causes the buckling modes. Usually, for a buckling problem, one often only 
observes the first mode (N = 1; the two transformation fronts are parallel). 
However, in the present problem, the second eigenvalue (N = 2) is very close to 
the first eigenvalue. Thus, if there is a slight disturbance, the second mode may 
appear. This may explain why the front can switch to an opposite inclination 
direction. 

We also point out that the solutions obtained above are also valid for a semi- 
infinite strip with the boundary conditions: 

V = 0, V x = 0, at X = +00, (6.30) 

if we restrict the spatial interval to < X < +00. Then, the right-half part of 
Figure 3 describes the different modes of a single transformation front, which 
is initially located near the left end. This corresponds to one experimental 
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situation in Shaw & Kyriakides (1998) (see Fig. 1 of that paper). 



7 Conclusions 

In order to better understand several instability phenomena observed in ex- 
periments in a thin SMA strip during the process of stress-induced phase 
transitions, we carry out an analytical study. Mathematically, it is a very chal- 
lenging problem as one needs to study the solution bifurcations of nonlinear 
partial differential equations in order to capture the instability phenomena. We 
start from the formulation of the three-dimensional field equations. By using 
the smallness of the thickness and the maximum strain, through a methodol- 
ogy which combines series expansions and asymptotic expansions, we derive 
the two-dimensional asymptotic equations, which take into account the lat- 
eral deformation and satisfy the traction-free boundary conditions up to the 
right order. Then, by further using the smallness of the width, we derive two 
one-dimensional asymptotic normal form equations for the phase transition 
problem in a thin strip. These two equations are analyzed and we manage to 
obtain some interesting analytic solutions for an infinite long strip under free- 
end boundary conditions through the WKB method. Our analytical results 
capture several instability phenomena observed in experiments successfully. 
It is shown analytically that the inclination of the transformation front is 
a phenomenon of localization-induced buckling (or phase-transition-induced 
buckling as the strain localization appears due to the phase transition). Also, 
it is demonstrated that there exists a second mode with a stress eigenvalue 
very close to that the first one. Thus, a slight disturbance could cause a switch 
from the first mode to the second mode. This, in turn, implies a switch of the 
inclination of the transformation front to an opposite direction, which ex- 
plains the orientation instability. Our results also reveal more explicitly the 
important role played by the thickness of the strip and show that the width 
influences the instability phenomena through the thickness-width ratio rather 
its magnitude. In literature, whether the well-known phenomenon of the Lud- 
ers band in a mild steel is caused by microscopic effects or macroscopic effects 
is still not completely settled issue. Due to the similarities between stress- 
induced transformations in a SMA and the development of Luders band in 
a mild steel, the present results also provide a strong mathematical evidence 
that the formation of the Luders band is a phenomenon due to macroscopic 
effects. 

Finally, we point out that as a by-product of the present study we have also 
provided a mathematically consistent derivation of the classical Euler buckling 
equation, without using the ad hoc hypothesis that the bending moment is 
proportional to the curvature. 
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Appendix: Incremental elastic moduli 



For initially isotropic material, in the case that there are no prestresses, 
$ should be a function of the principle stretches Ai, A2 and A3, namely $ = 
$(Ai,A 2 ,A 3 ). Denote by = J^|a 1 =a 2 =a 3 =i > $1 = $2 = $3 should vanish 
since there are no prestresses. The non-zero first order incremental elastic 
moduli can be written as 



A i = a im = $n> 
A 2 = a\ 122 = $12, 

A 3 = a\ 212 = -(A 1 -A 2 ), 
A 4 = a\ 221 = A 3 , 

where A 2 and A 3 are exactly the Lame's constants for infinitesimal strain. 
There are only two independent constants among A^ i = 1, 2, 3, 4. 

The non-zero second order incremental elastic moduli can be written as 



B 1 


2 

~~ a iiiin 


= $111, 


B 2 


2 

— a 111122 


= $112, 


£3 


2 

— °112233 


= $123, 




2 

~~ a lll212 


= ^{2A 2 + 2A 3 + B 1 -B 2 


B 5 


2 

— a 331212 


= ±(A 2 + B 2 -B 3 ), 


B 6 


2 

— a 121323 


= -{B A - B 5 ), 


B 7 


2 

— "111221 


= B 4 -A 2 - A 3 , 


B 8 


2 

— a 331221 


= B 5 - A 2 , 


B 9 


2 

— a 123123 


= B 6 - A 3 . 



There are only three additional independent constants among Bi,i = 1 ~ 9. 
The non-zero third order incremental elastic moduli can be written as 



Ci 


3 

= "llllllll 


= $1111, 


c 2 


3 

= a 11111122 


= $1112, 


c 3 


3 

= a llll2222 


= $1122, 


C4 


3 

= a llll2233 


— $1123, 


c 5 


3 

= a 11111212 
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C7 = o? 1221212 = -^(6^2 + 6A 3 - QB 2 - d - 2d + 3d), 

d = a 11221313 = — 4(^2 ~ B 2 — B 3 — C 2 + d), 

= ^(6A 2 + 6A 3 + 6B1 - 65 2 + d - 4C 2 + 3d), 

o 



Cg — a 12121212 

~ 3 



3 1 

do = a 12121313 = 0^9! 



C11 — a iii2i323 _ ~2A^ J ^ 1 ~ ^ — C\ + C 2 + 3C3 — 3C4), 
C12 = a?mi22i = ^(^2 + 6A 3 - SB, - W 2 + 2d - 2d), 

Cl3 = a 11112332 — ~2^ 2 _ ^ 3 + ^ 4 )' 

C14 = a? 1221221 = ^(6A 2 + QA 3 - QB 2 + d + 2C 2 - 3d), 

Cl5 = a 11221331 = 4(^2 — -B2 — B 3 + C 2 — C 4 ), 

Cie = 0L121221 = "^(6^2 + 6A 3 - d + 4C 2 - 3C 3 ), 

C17 = ^2211221 = l( QA 2 + 6 A 3 - 2Si + 2£? 2 + d - 4C 2 + 3d), 

o 



Cl8 — a 12121331 — 0^16, 



C19 = a?22ii33i = ^(6^i - 3Sx - 35 2 + QB 3 + d - 4C 2 + 3d), 

do = O11123123 — 24 + ^ 2 — 3-Bx + 3-B 3 + d — d — 3d + 3d), 

dl = a 11123132 = C20 — ^(^1 + ^2 + -Bi — B 2 ), 

C22 = ^2122323 = 2^(6^2 - 6A 3 + 125 2 - 125 3 + d - 4d + 3d)- 

There are only four additional independent constants among Ci,i = 1 ~ 22. 
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